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Abstract. The authors consider a mathematical model for the coupled atmosphere- 
ocean system, namely, the coupled quasigeostrophic flow-energy balance model. 
This model consists of the large scale quasigeostrophic oceanic flow model 
and the transport equation for oceanic temperature, coupled with an atmo- 
spheric energy balance model. After reformulating this coupled model as a 
random dynamical system (cocycle property), it is shown that the coupled 
quasigcostrophic-energy balance fluid system has a random attractor, and un- 
der further conditions on the physical data and the covariance of the noise, the 
system is ergodic, namely, for any observable of the coupled atmosphere-ocean 
flows, its time average approximates the statistical ensemble average, as long 
as the time interval is sufficiently long. 



1. Mathematical model 

We consider large scale geophysical flows modeled by the quasigeostrophic flow 
equation in the horizontal xy— plane, in terms of vorticity q(x,y,t), and the trans- 
port equation for the oceanic temperature T[x, y, t), coupled with the atmospheric 
energy balance equation proposed by North and Cahalan ^B] for the air tempera- 
ture Q(x, y, t), on the domain D = {(x, y) : < x < I, < y < I}: 

9 t =A9 - (a + 9) + S a (x, y) - b{y){S {x, y) + 6 - T(x, y)) + w, 
dq 

(1) -g^ =vAq - rq + Pr Ra d y T - J(tp,q + (3y), 

f =AT-J(T^), 

where -0(ir,y,£) is the stream function, (3 > is the meridional gradient of the 
Coriolis parameter, v > is the viscous dissipation constant, and r > is the 
Ekman dissipation constant. Furthermore, 

q(x,y,t) = Aip(x,y,t) 

is the vorticity, a > is a constant parameterizing the effect of the earth's long- 
wave radiative cooling, b(y) is the latitudinal fraction of the earth covered by the 
ocean basin, Pr is the Prandtl number, and Ra is the Rayleigh number. Note 
that S a (x,y) and S a (x,y) are empirical functions representing the effects (on at- 
mosphere and ocean, respectively) of the shortwave solar radiation. Moreover, 
J(<7, h) — g x h y — g y h x is the Jacobian operator and A = d xx + d vv is the Laplacian 
operator. All these equations are in non-dimcnsionalizcd forms. The fluctuating 
noise w(x,y,t) is usually of a shorter time scale than the response time scale of 
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the air mean temperature. So we neglect the autocorrelation time of this fluc- 
tuating process and thus assume that the noise is white in time. The spatially 
correlated white-in-time noise w(x,y,t) is described as the generalized time deriv- 
ative of a Wiener process w(x,y,t) defined in a probability space (f2,F, P), with 
mean zero and covariance operator Q. For geophysical background of similar cou- 
pled atmosphere-ocean models, see [El El- Recently, a few authors have considered 
the randomly forced quasigeostrophic equation, in order to incorporate the impact 
of uncertain geophysical forces ( 0)- 

The fluid boundary condition is no normal flow and free-slip on the whole boundary: 



(2) if, = 0, q = 0. 

The flux boundary conditions are assumed for the ocean temperature T and air 
temperature 0: 

5e ar n 

3 7T = tt = °> 

On On 

with n is the outer unit normal vector on boundary. By equation and the 
boundary condition, we find that J D T(x, y)dxdy — constant, and thus we assume 
that J D T(x, y)dxdy — without loss of generality. 

This is a coupled system of both deterministic and stochastic partial differential 
equations. 

This paper is organized as follows. In the next section, we introduce basic concepts 
in random dynamical systems. Then we reformulate the coupled fluid system 0) as 
a random dynamical system with help of a Ornstein-Uhlenbeck random stationary 
process in §3. After obtaining basic estimates for the system in §4, we show in 
§5 that the coupled fluid system admits a random attractor. Under further 
conditions on the physical data and the covariance of the noise, we show that the 
system is actually ergodic, namely, for any observable of the coupled atmosphere- 
ocean flows, its time average approximates the statistical ensemble average, as long 
as the time interval is sufficiently long. 

2. Random dynamical systems 

In order to investigate the long time dynamics of the coupled fluid system (Q under 
the influence of random forces, we need some appropriate concepts and tools from 
the theory of random dynamical systems. 

A random dynamical system consists of two components. The first component is a 
driven flow (O, P, 9) as a model for noise, where (fi, J 7 , P) is a probability space 
and 9 is a T <E> B(M.),J- measurable flow: we have 



9q — id, 9t+ T — 9t o 9 r =: 9t9 T 

for t, t € K. To express that the noise is stationary and "chaotic" , the measure P 
is supposed to be ergodic with respect to 9. The second component of a random 
dynamical system is a £?(R + ) <g> T <£> B(H), £>(iT)-measurable mapping ip satisfying 
the cocycle property 

<p(t + r, u, x) = ip(t, 9 t uj, ip(r 7 w,x)), ip(0, u, x) = x, 

where the phase space H is a separable metric space and x is chosen arbitrarily in 
H. We will denote this random dynamical system by symbol (p. 
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A standard model for a spatially correlated white-in-time noise is the generalized 
time derivative of a two-sided Brownian motion or Wiener process w(x,y,t). Let 
U be a separable Hilbert space with scalar product (•, •) and the induced norm || • ||. 
For stochastic partial differential equations containing such a noise, an appropriate 
or canonical probability space is 

(C (R, U),B(C (R, £/)), P), 

where Co(M, U) is the space of continuous functions on R, which take zero value at 
time zero. This space is given the compact-open topology (i.e., uniform convergence 
on compact intervals in R). This topology is metrizible as it can be generated by 
the complete metric 

^ 1 d n (g 1 ,g 2 ) 



d{gi,92) 



^ 2 n l + d n ( 9l ,g 2 y 



where d n (g 1 ,g 2 ) = max| t |<„ \\gi{t) - g2{t)\\, for gi,g 2 in C (R, U). Thus we have 
open balls or open sets in Co(R 7 U), and B(Co(R, U)) is the corresponding Borel 
cr-algebra. Suppose the Wiener process w has covariance operator Q on U . Let P 
denote the Wiener measure with respect to Q. Note that P is ergodic with respect 
to the Wiener shift 9 t : 

(4) 6 t u = u{- + t)-u{t), for w G C (R,U). 

A major example of a random dynamical system is a random differential equation. 
For example, let us consider the following evolution equation in some Hilbert space 



(5) — = f(u,0 t u), u(Q)=x, 

at 

over some metric dynamical system (f2, T , P, 0). If JSJ is well-posed for every u> G fi 
and solutions u(t, w,x) depend measurably on (t,u>,x), then the operator 

tp : (t, uj, x) — y u(t, ui; x) 

defines a random dynamical system (cocycle) ip. For detailed presentation of ran- 
dom dynamical systems we refer to the monograph by Arnold 1 . 

Motivated by deterministic dynamical systems we introduce several useful concepts 
from the theory of random dynamical systems. 

A closed set B(cu), depending on u, in a separable Hilbert space H is called random 
if the distance mapping lu — > sup^g^^ ||x— u\\h is a random variable for any y E H. 
A random dynamical system is called dissipative if there exists a random set B that 
is bounded for any lo and that is absorbing: for any random variable x(lo) G H 
there exists a t x (uj) > such that if t > t x (oj), then 

(p(t,oj,x(uj)) G B(6 t uj). 

In the deterministic case ((p is independent of w) the last relation coincides with 
the definition of an absorbing set. In the case of partial differential equations of 
parabolic type, due to the smoothing property, it is usually possible to prove that a 
dissipative system possesses compact invariant absorbing sets. For more details, see 
Temam [25]. page 22f. Hence for a system of stochastically forced parabolic partial 
differential equations, such as the stochastic two-layer fluid system introduced in the 
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last section, we usually consider the random set B(w) to be compact. In addition, 
we will assume that B(u>) is forward invariant: 

(p(t,u,B(u>)) C B(d t u), t > 0. 

In the following we also need a concept of tempered random variables. A random 
variable x is called tempered if 

t -» \x(9 t uj)\ 

is subexponentially growing: 

limsup — — flX^I! — o a.s. 

t^±oo \t\ 

where 

log + {x) = max{0,log(x)} 

This technical condition is not a very strong restriction because the only alternative 
is that the above limsup is oo, which describes the degenerate case of stationarity; 
see Arnold P, page 164 f. 

3. COCYCLE PROPERTY 

In this section we will show that the coupled fluid model 0) defines a random 
dynamical system (cocycle property). The cocycle property essentially comes from 
the well-posedness. 

Now we re- formulate the model such that appropriate tools of the theory of random 
dynamical system can be applied to analyse the coupled atmosphere-ocean model 
under a random wind forcing. For the following we need some tools from the theory 
of partial differential equations. 

Let H 1 (D) be the Sobolev space of functions on D with first generalized derivative 
in Z/2(-D), the function space of square integrable functions on D with norm and 
inner product 

\\u\\l 2 =(J \u(x)\ 2 dL)\ , (u,v)l 2 = J u(x)v(x)dD, u, v € L 2 (D). 

The space H 1 (D) is equipped with the norm 

IMIff 1 = \W\\l 2 + \\d y u\\ L2 + \\d z u\\ L , 2 . 

o 

Motivated by the zero-boundary conditions of q we also introduce the space H (D) 
which contains roughly speaking functions which are zero on the boundary dD of 
D. This space can be equipped with the norm 

(6) \\uU HD) = \\d y u\\ L2 + \\d z u\\L 2 . 

Another Sobolev space is given by H l (D) which is a subspace of H 1 (D) consisting 
of functions u such that J D udxdy — 0. A norm equivalent to the iJ^-norm on 
ij 1 (D) is given by the right hand side of ©■ For functions in L 2 (D) having this 
property we will write L2(D). 

For convenience, we introduce the vector notation for unknown geophysical quan- 
tities u = (0,q,T) 

We now take the linear differential operator from JJJ 
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(l + b(x,y))e 



Recall that the function < b(y) < 1. A is defined on functions that are sufficiently 
smooth. We also have the following boundary conditions from §.2 

d& _ dT 
dn dn 

ip\dD = 0, q\dD = 0. 
We introduce the phase space for our geophysical quantities 

H = L 2 (D) x L 2 (D) x L 2 {D) 

V = H X {D) x H^D) x H^D). 

After this preparation, we are able to write our problem as a stochastic evolution 
equation. 

Let w be a noise on L 2 {D) with finite energy given by the covariance operator Q 
of the Wiener process w(t) which is defined on a probability space P). For 

the vector 

W = (w,0,0), 

we rewrite the coupled atmosphere-ocean system as a stochastic evolution equa- 
tion on V': 

du 

(7) — + Au = F{u) + W, u(0) — uq G H, 

at 

where w is a white noise as the generalized temporal derivative of a Wiener process 
w with continuous trajectories on K and with values in L 2 (D). And F(u) is defined 
as the follow: 

/ -a + S a (x, y) - b{y)S (x, y) + b(y)T(x, y) 
F(u) = -rq + PrRad y T - J(i/j, q + (3y) 
V -J(T,V) 

Sufficient for this regularity is that the trace of the covariance is finite with respect 
to the space L 2 (D): ti'L 2 Q < oo. In particular, we can choose the canonical prob- 
ability space where the set of elementary events O consists of the paths of w and 
the probability measure P is the Wiener measure with respect to covariance Q. 

In the following, we need a stationary Ornstein-Uhlenbeck process solving the linear 
stochastic equation on D 

/ox dz 
8 — -+A lZ = w 

at 

where A\ = —A + (1 + b(y)) is the linear operator with zero Neumann boundary 
condition at dD and with zero initial condition. 

Lemma 3.1. Suppose that the covariance Q has a finite trace : tiL 2 Q < oo. Then 
0) has a unique stationary solution generated by 
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Moreover, Z(u>) = (z(oj), 0,0) is a random variable in V. 

For the proof we refer to Da Prato and Zabczyk ^7], Chapter 5, or Chueshov and 
Scheutzow @|. 

For our calculations it will be appropriate to transform (J7J) into a differential equa- 
tion without white noise but with random coefficients. We set 

(9) v:=u-Z 
Thus we obtain a random differential equation in V 

dv 

(10) — + Av = F(v + Ztftu)), v(0) = v EH. 
dt 

Equivalently, we can formulate the equation l|10l) using test functions 

f t (v(t),Q +a(v(t),Q = (F(v{t)+Z{0tu)),Q for all £ G V. 

We have obtained a differential equation without white noise but with random 
coefficients. Such a differential equation can be treated sample-wise for any sample 
oj. Hence it is simpler to consider (|10|) than to study the stochastic differential 
equation J7J directly. We are looking for solutions in 

veC([0,r};H)nL 2 (0,r; V), 

for all r > 0. If we can solve this equation then u := v + Z defines a solution version 
of 0. For the well posedness of the problem we now have the following result. 

Theorem 3.2. fWell-PosednessJ For any time r > 0, there exists a unique 
solution of $10\) in C([0, r]; H) n £2(0, r; V). In particular, the solution mapping 

R+ x D, x H 3 (t,oj,v ) -> v(t) e H 

is measurable in its arguments and the solution mapping H 9 vo — > v{t) G H is 
continuous. 



Proof. By the properties of A and F the random differential equation (|10|l is essen- 
tially similar to the 2 dimensional Navier Stokes equation. Hence we have existence 
and uniqueness and the above regularity assertions. □ 

On account of the transformation © , we find that (7J also has a unique solution. 
Since the solution mapping 

IR+ x O x H 9 (f, u, vo) — > v(t, to, v ) =: tp(t, to, v ) G H 

is well defined, we can introduce a random dynamical system. On f2 we can define 
a shift operator t on the paths of the Wiener process that pushes our noise: 

w(-,6 t w) = w(- +t,uj) — w(t,uj) for t G R 

which is called the Wiener shift. Then {0 t }tes. forms a flow which is ergodic for the 
probability measure P. The properties of the solution mapping cause the following 
relations 

ip(t + t, to, u) = ip(t, 9 t lu, (p(r, oj, u)) for t, r > 
ip(0, oj, u) = u 
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for any to S and u G H. This property is called the cocycle property of (p which 
is important to study the dynamics of random systems. It is a generalization of 
the semigroup property. The cocycle if together with the flow 6 forms a random 
dynamical system. 



4. DlSSIPATIVITY 

In this section we are going to show that the coupled atmosphere-ocean system Q 
is dissipative, in the sense that it has an absorbing (random) set. This definition 
has been used for deterministic systems [23]. This means that the solution vector 
v is contained in a particular region of the phase space H after a sufficiently long 
time. Dissipativity will be very important for understanding the asymptotic dy- 
namics of the system. This dissipativity will give us estimate of the atmospheric 
temperature evolution under oceanic feedback. Dynamical properties that follow 
from this dissipativity will be considered in the next section. In particular, we will 
show that the coupled atmosphere-ocean system has a random attractor, has finite 
degree of freedom, and is ergodic under suitable conditions. 

Let 

(11) v = (Q,T), Q = e-z. 
By integration by parts and Poincare inequality, we have 

(12) M 2 . > 

-^Wh frill, + 6a 2 + 6US.Hk + 6||S ||| 2 . 

Here Poincare inequality Ao||w|| 2 < ||Vm|| 2 (Ao > 0) is used. Here and in the follow- 
ing we stress that < b(y) < 1. 
For q, we have 

(13) ^- t \\l\\h = -"\\Vq\\h - r\\q\\h + PrRa J^d y Tqdxdy - JJ^qdxdy. 
Under the following condition 

B 2 \l\ 2 

(14) Avr > PJlL 

followed the argument of |12| . we have for some positive constant a > 

d f 
^ t h\\h < -a\\q\\ 2 L2 + PrRa j^d y Tqdxdy, 

by the Cauchy-Schwarz inequality, we get 

d , „, . „ „, Pr 2 Ra 2 
dt a 
Collecting all these estimates, we have for some positive constants ot\ and C\ 



(15) -rMWh < -<*\\Q\\h + — llVTHi, 



|^l|B|lk + l|T||k + ^|| ? ||I 2 ) + « 2 ||V,f 

<-«i(^lie||^ + ||T||| 2 + ^^||# i2 ) + c 1 
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and 



(17) fwvvirnidT^cwvmi 



+ C 5 (t - s), for any < s < t. 



s 



We now get the dissipativity of v. Roughly speaking dissipativity means that all 
trajectories of the system move to a bounded set in the phase space. For a random 
system we have the following version of dissipativity. 

Definition 4.1. A random set B = {BIlu)}^^ consisting of closed bounded sets 
B(u) is called absorbing for a random dynamical system (p if we have for any 
random set D = {_D(w)} we o, D(ui) £ H bounded, such that t — > sup ye£) ( etaj ) \\y\\ii 
has a subexponential growth for t — > ±oo 

ip{t,ui,D{u}))C.B{9 t u}) fort>tv{D,uj) 
l ' M ^(M- t w,DM)cBH fort> t (D,u). 

B is called forward invariant if 

ip(t, u, u ) G B(9 t uj) if u € B(uj) for t > 0. 

From l|16|l . we could get the existence of absorbing set B(lo) = {v 6 H, || v j| ^ < 
2C3}. Suppose that t — » |jwo(^-t^)||lf growths not faster than subexponential. 
Then we have the subexponential growth for B(lo), and we get 

Lemma 4.2. The random set B(uj) = B(0,R(u)) = {v £ H, \\v\\ 2 H <R = 2C 3 } is 
an absorbing and forward invariant set for the random dynamical system generated 

by m>. 

For the applications in the next section we need that the elements which are con- 
tained in the absorbing set satisfy a particular regularity. To this end we introduce 
the function space 



where s£l. The operator A s is the s-th power of the positive and symmetric op- 
erator A. Note that these spaces are embedded in the Slobodeckij spaces H s , s > 0. 
The norm of these spaces is denoted by || • \\h s - This norm can be found in Egorov 
and Shubin ^Hl, Page 118. But we do not need this norm explicitly. We only 
mention that on H s the norm || • || s of H s is equivalent to the norm of TL S for < s, 
see [T5] . 

Our goal is it to show that v(1,uj, D) is a bounded set in H s for some s > 0. This 
property causes the complete continuity of the mapping v(l,u>, •). We now derive 
a differential inequality for By the chain rule we have 



H s :={ueH: \\uf s := \\A*uf H < 00} 



i(t||«(*)[[2) = IK*)H2 + 4lK*)H2- 



Note that for the embedding constant ce, s between 7i s and V 




for s < 1 
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such that the left hand side is bounded if the initial conditions vq are contained in 
a bounded set in H . The second term in the above formula can be expressed as 
followed: 

tj t {^v, Aiv) H = 2t(^w, A s v) H = -2t(Av, A s v) H + 2t(F(v + Z(6 t w)), A s v) H - 
We have 

(Av,A°v) H = \\A^v\\ H = \\v\\ 2 1+s . 
For terms which including Jacobian operator, we apply some embedding theorems, 
see Temam [221 Page 12, then taking F3 as a example we have for a constant C5 > 

{F 3 (v),Oh < C 6 |k||m 1 + l||^||m 2 + l||CI|m3, C £ #m 3 

where mi + ni2 + mi?, > 1 and < mj < 1. Here we use that D is of dimension 2. 
We then have for mi — 0, 7712 = s < 1 and m? = 1 — s 

|(F3(v) !J 4 s T) ff |< C7 ||T|| y ||^|| 1+s ||T|| 1+s . 

||V>(*)l|i+s is bounded by ci||g(t)||/f by the definition of ip(t) and |M*)IUoo(o,T;ir) < 
00. Hence we have for any e > a constant c$(e): 

\\T(t)r v + e\\T(t)\\l +s 

<C8(e)ll«llL(o,r;H)ll«(*)llv + e|l«(*)lli+.. 
where e is chosen sufficiently small. The other terms could be estimated similarly, 
we omit it here. Using (|T7jl . we could obtain \\v(t, u, Uo)|| s , < s < \ is bounded 
for t < T < 00, t > if v is contained in a bounded set. This allows us to write 
down the main assertion with respect to the dissipativity of this section. 

Theorem 4.3. For the random dynamical system generated by there exists a 

compact random set B — {_B(w)} w6 n which satisfies Definition ^. 1\ 

We define 

(19) B(ui) = <p(l,0_ia;,£(O,.R(0_iw))) C U s , < s < J. 

In particular, TL S is compactly embedded in H. 

5. Random dynamics 

In this section we analyse the dynamical behavior of the coupled atmosphere-ocean 
system However, it will be enough to analyse the transformed random dynam- 
ical system generated by {7J). By the transformation © we can take over all these 
qualitative properties to the system JTJl. 

We will consider following dynamical behavior: random attractors, atmospheric 
temperature evolution under oceanic feedback, and ergodicity. 
We first consider random climatic attractors. We recall the following basic concept; 
see, for instance, Flandoli and Schmalfufi 

Definition 5.1. Let if be a random dynamical dynamical system. A random set 
A = {A^uj)}^^ consisting of compact nonempty sets A(uj) is called random global 
attractor if for any random bounded set D we have for the limit in probability 



(P) lim dist H (v(t,uj,D(u)),A{0 t Lj)) = 

t— »oo 
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and 

tp(t,uj,A(u)) =A{0 t u) 

any t > and u> G Q. 

The essential long-time behavior of a random system is captured by a random 
attractor. In the last section we showed that the dynamical system ip generated by 
(01 is dissipative which means that there exists a random set B satisfying . In 
addition, this set is compact. We now recall and adapt the following theorem from 

EH 

Theorem 5.2. Let ip be a random dynamical dynamical system on the state space H 
which is a separable Banach space such that x — > tp(t, U),x) is continuous. Suppose 
that B is a set ensuring the dissipativity given in definition ^. 1\ In addition, B has 
a subexponential growth (see Definition and is regular (compact). Then the 
dynamical system ip has a random attractor. 

This theorem can be applied to our random dynamical system ip generated by the 
stochastic differential equation (0. Indeed, all the assumptions are satisfied. The 
set B is defined in Theorem 14.31 Its subexponential growth follows from B(ui) C 
B(0, R(ui)) where the radius R(u>) has been introduced in the last section. Note 
that ip is a continuous random dynamical system; see Theorem 13.21 Thus ip has 
a random attractor. By the transformation this is also true for the original 
coupled atmosphere-ocean system. 

Corollary 5.3. (Ttandom Attractor ) The coupled atmosphere- ocean system 0) 
has a random attractor. 

Now we consider random fixed point and ergodicity. We can do a small modification 
of QJ. This modification is given when we replace AT by is AT where v > is 
viscosity. Under particular assumptions about physical data in (IJ , we can show 
that the behavior of our dynamical system is laminar. For a stochastic system, 
this means that after a relatively short time, all trajectories starting from different 
initial states show almost the same dynamical behavior. This can be seen easily if 
a, S a and S are zero, there is no noise and v is large. We will show that a laminar 
behavior also appears when a, S a and S are small in some sense. 
Mathematically speaking, laminar behavior means that a random dynamical system 
has a unique exponentially attracting random fixed point. 

Definition 5.4. A random variable v* : Q — ► H is defined to be a random fixed 
point for a random dynamical system if 

ip(t,U,V*(u))=V*(e t L0) 

for t > and uj G tt. A random fixed point v* is called exponentially attracting if 

lim \\ip(t,u,x)-v*(6 t Lu)\\ H = 

t — >oo 

for any x £ H and 10 £ fl. 

Sufficient conditions for the existence of random fixed points are given in Schmalfufi 
We here formulate a simpler version of this theorem and it is appropriate for 
our system here. 
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Theorem 5.5. (Ptandom Fixed Point Theorem,) Let ip be a random dynamical 
system and suppose that B is a forward invariant complete set. In addition, B has 
a subexponential growth, see Definition ^. 1\ Suppose that the following contraction 
conditions holds: 

(20) sup M^i) ~ <p(^)\\h < k{u) 

Vi =£v 2 £B(lj) \\Vi-V2\\h 

where the expectation of log k denoted by Elogfc < 0. Then ip has a unique random 
fixed point in B which is exponentially attracting. 

This theorem can be considered as a random version of the Banach fixed point 
theorem. The contraction condition is formulated in the mean for the right hand 
side of (|5U|) . 



Theorem 5.6. Assume that the physical data \a\, \\S a \\L 2 > \\S \\l 2 and the trace 
of the covariance for the noise tr#Q are sufficiently small, and that the viscosity 
v is sufficiently large. Then the random dynamical system generated by has a 
unique random fixed point in B. 

Here we only give a short sketch of the proof. Let us suppose for a while that B is 
given by the ball B(0, R) introduced in Lemma f4. 21 Suppose that the data in the 
assumption of the lemma are small and v is large. Then it follows that Ei? is also 
small. To calculate the contraction condition we have to calculate ||<p(l, w, v\{lo)) — 
ip(l, u>, 1)2(1-0)) \\ 2 H for arbitrary random variables v%, v 2 £ B. By the property of 
Jacobian operator we have that 

(J(gi,^i) - J(qi, ipi),qi - q2) < c 9 \\qi - g 2 ||^i + cio Iki II ^1 (jD) ||<?i - <?2||i 2 

where the constant C10 can be chosen sufficiently small if v is large. By direct 
estimate, we get 

^\\ip(t,uj,vi{uj)) - ip(t,u),v 2 (u)))\\ 2 H 

< (-a' + c w \\<p(t,(J, v 2 (uj)) \\v)\\ ip(t,u),vi(u))) - cp(t,u,v 2 {u)))\\% 

for some positive a' depending on v. From this inequality and the Gronwall lemma 
it follows that the contraction condition l|20() is satisfied if 

E sup cio / \\ip(t,uj,V2)\\vdt < a . 

u 2 G-B(w) JO 

But by the energy inequality this property is satisfied if the Ei? and E||z||y is suf- 
ficiently small which follows from the assumptions. 

Let now B be the random set defined in l|19() . Since the set B introduced in l|19() 
is absorbing any state the fixed point v* is contained in this B. In addition v* 
attracts any state from H and not only states from B. 

Corollary 5.7. (Unique Random Fixed Point,) Assume that the physical data 
\a\, \\S a \\L 2 ! II'S'oIIl 2 an d the trace of the covariance for the noise trnQ are suf- 
ficiently small, and that the viscosity v is sufficiently large. Then, through the 
transformation the original system has a unique exponentially attracting 
random fixed point u*(u>) = v*(ui) + Z(uS), where u = (<d,q,T). 



12 



AIJUN DU, JINQIAO DUAN, HONGJUN GAO, AND TAMAY OZGOKMEN 



The uniqueness of this random fixed point implies ergodicity. We will comment on 
this issue at the end of this section. 

By the well-posedness Theorem 13.21 we know that the stochastic differential equa- 
tion J7J for the coupled atmosphere-ocean system has a unique solution. The solu- 
tion is a Markov process. We can define the associated Markov operators T(t) for 
t > 0, as discussed in [2D [2D] • Moreover, {T(i)} t > forms a semigroup. 
Let M 2 be the set of probability distributions fi with finite energy i.e., 



\u\\ 2 H dfi(u) < 



H 



Then the distribution of the solution u(t) (at time t) of the stochastic differential 
equation J7J is given by 

2"(*)Mo, 

where the distribution of the initial data is contained in M 2 . 

We note that the expectation of the solution ||it(i)||jj can be expressed in terms of 

this distribution T(t)fio: 

nHt)\\ 2 H = f \\ u \\ 2 H dT{t)^. 

J H 

We can derive the following energy inequality in the mean, using our earlier esti- 
mates in ifTB)) and lfT7| : 

Theorem 5.8. The dynamical quantity u = (O, q, T) of the coupled atmospheric- 
ocean system( 1 ) satisfy the estimate 

¥]\u(t)f H +cM f ||tt(T)||^dr<E||u ||^ + *cii+*tr£ a Q J 
Jo 

where the positive constants Cn and a depend on physical data a, ||5' a ||i2, ||S' ||l2, 
Pr and Ra. 

By the Gronwall inequality, we further obtain the following result about the as- 
ymptotic mean-square estimate. 

Corollary 5.9. (Asymptotic Mean-Square Estimate ) For the expectation of 
the dynamical quantity u — (0,g, T) of the coupled atmospheric- ocean system (1), 
we have the asymptotic estimate 

en + tr L2 Q 



limsupE||u(i)||^ = limsup / \\u\\ 2 H dT (t) fi < 

t-^oo t — >oo J H 



Cl2 

if the initial distribution fj,Q of the random initial condition Uo(u>) is contained in 
M 2 . Here c\2 > also depends on physical data. In particular, we have asymp- 
totic mean-square estimate for the atmospheric temperature evolution under oceanic 
feedback 

(21) limsupEHeil^ < Cll+trL2Q . 

t— >oo C12 

Thus the atmospheric temperature Q(y,t), as modeled by the coupled atmosphere- 
ocean system 0), is bounded asymptotically in mean-square norm in terms of phys- 
ical quantities such as the trace of the covariance operator of the external noise, the 
earth's longwave radiative cooling coefficient a, and the empirical functions \\S a \\L 2 
and || S 1|^2 representing the latitudinal dependence of the shortwave solar radiation, 
as well as the Prandtl number Pr and the Rayleigh number Ra for oceanic fluids. 
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By the estimates of Theorem 15.81 we are able to use the well known Krylov- 
Bogolyubov procedure to conclude the existence of invariant measures of the Markov 
semigroup. 

Corollary 5.10. The semigroup of Markov operators {T(t)}t>o possesses an in- 
variant distribution /ij in M 2 : 

T{t)ni = fa for t>0. 

In fact, the limit points of 

1 1 Jo ) t>0 

for t — > oo are invariant distributions. The existence of such limit points follows 
from the estimate in Theorem 15. 81 

In some situations, the invariant measure may be unique. For example, the unique 
random fixed point in Corollary 15.71 is defined by a random variable u*(lu) — 
v*(lo) + z(u>). This random variable corresponds to a unique invariant measure 
of the Markov semigroup. More specifically, this unique invariant measure is the 
expectation of the Dirac measure with the random variable as the random mass 
point 

Because the uniqueness of invariant measure implies ergodicity |18| . we conclude 
that the coupled atmosphere-ocean model Q is ergodic under the suitable condi- 
tions in Corollarv l5 . 71 for physical data and random noise. We reformulate Corollary 
15.71 as the following ergodicity principle. 

Theorem 5.11. (T^rgodicityj Assume that the physical data \a\, \\S a \\^2, ||S , ||^2 
and the trace of the covariance for the noise tr hQ are sufficiently small, and that 
the viscosity v is sufficiently large. Then the coupled atmosphere- ocean system 0) is 
ergodic, namely, for any observable of the coupled atmosphere- ocean flows, its time 
average approximates the statistical ensemble average, as long as the time interval 
is sufficiently long. 
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